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The exchange contribution, ni(k, 0), to the static dielectric function in the electron liquid is eval¬ 
uated exactly. Expression for it is derived analytically in terms of one quadrature. The expression, 
as presented in Eq. (3) in the Introduction, turns out to be very simple. A fully explicit expression 
(with no more integral in it) for ni(k, 0) is further developed in terms of series. Equation (3) is 
proved to be equal to the expression obtained before under some mathematical assumption by Engel 
and Vosko, thus in the meanwhile putting the latter on a rigorous basis. The expansions of ni(k, 0) 
at the wavectors of A: = 0, k = 2kp, and at limiting large k are derived. The results all verify those 
obtained by Engel and Vosko. 

PACS numbers: 71.15.Mb, 71.10.-w, 71.45.Gm 


I. INTRODUCTION WITH CONCLUDING 
REMARKS 

The dielectric function in the homogeneous electron 
liquid is determined by Hi 


e(k, w) = 1 — i;(A:)n(k, w); (1) 


n(k, w) is the proper (irreducible) linear response func¬ 
tion. [v{k) is the Fourier transform of the Coulomb po¬ 
tential.] Many electronic properties in solids and liquids 
can be described from the dielectric function. Rapid de¬ 
velopments in the techniques of inelastic x-ray scattering 
and electron energy-loss spectroscopies in recent years 
now enable us to measure these dielectric function 
related properties with unmecedented accuracy. 

Theoretically, Lindhard [1| obtained, about six decades 
ago, the explicit expression for no(k, w), the zeroth-order 
of n(k, w) in terms of Coulomb interaction, in the homo¬ 
geneous electron gas; wherefore it is also known as Lind¬ 
hard function. Keeping no(k,w) only in IM. (1) yields 
the random phase approximation (RPA) HIH for e(k, w). 
An enormous amount of effort had been devoted to the 
study of the correction beyond RPA since the pioneer¬ 
ing proposals made by Hubbard (and DuBois [HI). 
The correction beyond RPA is fully due to the exchange- 
correlation effects, which is also termed, according to 
Hubbard, as local field correction G(k, w): 


G(k, ui) = v{k) 


-1 


_n(k,w) no(k,w) 


( 2 ) 


Needless to say, the lowest order correction, ni(k, w), 
to the Lindhard function, which is traditionally also 
known as the exchange contribution to n(k, w), had 
been of major research interest [l2l - E^ . In some of 
the researches, not strictly ni(k, w) but variants of it 
were studied as well. Higher order contributions beyond 
nn(k, a;)-|-ni (k cc) had also been studied, for instance, 
in Refs. . About two and a half decades ago, En¬ 

gel and Vosko obtained an analytical expression [Eq. 
(29) in Ref. [^] for ni(k, 0) in terms of one quadra¬ 
ture. That was definite progress, and one must be aware 
that the original form for ni(k, 0) is in six-fold integral 
[see Eq. ([5]) in the next section]. Indeed, most of the 
work we mentioned above relied on numerical procedures 
in one way or another. Unfortunately the derivation by 
Engel and Vosko can not be fully regarded as rigorous in 
that it critically relies on some mathematical assumption. 
But convergence to the numerical results for Hi (k, 0) 
and many right analytical features all strongly suggest 
the correctness of their expression. Very encouragingly 
in this connection, Glasser [2^ confirmed Eq. (29) of 
Ref. 23| with the aid of analytical computer programs. 
Several of the leading terms in the small wavevector ex¬ 
pansion of ni(k, 0) obtained from it were also confirmed 
later analytically by Svendsen and von Barth [^ . 


In this paper we report an exact evaluation of Hi (k, 0). 
The main result of the present work is the following ex¬ 
pression for ni(k, 0), 


ni(k,0) = 


2 2 
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with a = 1 — fc/2, 6 = 1-1- k/2, and k in units of kp- 
The result for ni(k, 0) can be claimed to be surprisingly 


simple. Equation ([3]) is proved to be equal to the one 
given by Engel and Vosko 23| mentioned above; but in 
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the present work it is derived rigorously. The expansions 
of it at fc = 0, A: = 2, and at limiting large k are also de¬ 
rived; they fully confirm the corresponding ones obtained 
in Ref. [ 23 . 

The remaining integral in Eq. © is a straightforward 
numerical task. But we have an interesting alternative 
which is to further carry out analytically the remaining 
integral in terms of series. In this way, a fully explicit 
expression for ni(k, 0) is developed, which is shown in 
Eq. (IMl) in Sec. V, with the quantity I(k) being related 
to ni(k, 0) as in Eq. ([H]), and C(3) being Riemann’s zeta 
function. This series representation should be suitable 
for both of numerical and analytical applications. 

The details of our derivation are given in Sections II- 


IV; supporting materials are relegated to Appendices A 
and B. The final result is given in Sec. V in which the 
above-mentioned series representation is also developed. 
The proof of the equivalence of our result and Eq. (29) in 
Ref. [ 2 ^ is given in Sec. VI. The expansions of IIi (k, 0) 
in limiting cases are derived in Sec. VII. 


II. REDUCTION OF ni(k,0) TO 
TWO-DIMENSIONAL INTEGRAL 

Explicit expression for Hi (k, ca) can be obtained with 
the diagrammatic techniques [ll| (see also Refs. [3 Hi): 


J 


ni(k,w) =2 


dp dp p^) ^p+k)(^p^ ^p^+k) 


(27r)3 (27r)3 


OJ LUp — Wp+k T *0“*" 


UJ -f UJp — Wp+k T *0^ ^ + <^p' — Wp'+k T *0^ 


(4) 

where Wp = /2m and rip is the Fermi-Dirac distribution function. (We put /i = 1 in this paper.) At a; = 0, one gets 


2 2 f f ^P' 

ni(k,0) = - 8 ^m e J J ^^np_k/2np._k/2 


p' k 


p • k p'• k 


IP-P' 


.(5) 


The Hi (k, 0) in Eq. ([5]) is equal to the corresponding one 
defined in Ref. [231 ■ We also define / {k) as 

/(fc) = ^ni(k,o). (6) 

We note that ni(k, 0) depends only on the magnitude k 
in a uniform system. /(fc) so defined has the well-known 
property of /(fc = 0) = — 1 (see Sec. VII). 

We shall perform our calculation in the cylindrical 
coordinates. The advantage of the cylindrical coordi¬ 
nates in calculating Hi (k, uj) had been appreciated be¬ 
fore [3, Hll, [ 2 ^ 13.1^ . The integrals over the azimuthal 
angular variables of p and p' can be readily carried out. 
After that one obtains 


can be further carried through. The resulting expression 
for L{u) is 

L{u) = -[Li(m) — (A-I-A')(l-I-2 In |2u|) — n] 

-j- A In \Li (u) -t- A — A -t- u| 

-I- A'In |Li(u)-I-A — A'-|-m|, (9) 

where 

Li{u) = \/{X' — X + u)^ + AXu 

= + {X- X'Y + 2(A -f X')u. (10) 

Substituting Eq. (|H]) into Eq. one obtains 


J P^L{I3% 

(7) 

where a = z + z', /3 = z — z'. The same notations as 
those in Ref. [ 2 ^ are adopted. The function L in Eq. 
© is defined as 


L{u) 



dx' 


1 

\/{x — x'Y + 2{x + x')u ’ 


( 8 ) 

where X = {a + z){b — z) and A' = (a -I- z'){b — z'). The 
two-dimensional integral on the right hand side of Eq. © 


m = 


16fc2 




2=0 


( 11 ) 


where 


Jo = -4 / [ dzdz'^[{X + X'){2\n2 + l)+a^ + P^]; 

J —a J 

( 12 ) 

Ji= [ [ dzdz'^[a^L,{a‘^) - (13) 

J-aJ ^ ^ 
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J2 =2 


dzdz'- 


{ A[a2 In \Li{a'^) + A' - A + a^| - In \Li{P^) + A' - A + 


+ A'K In \Li{a^) + X - X' + a^\- \n + X - X'+ (3^\]}, 


(14) 


or, by the use of the fact that both of Li(a^) and Li(/3^) remain unchanged [in virtue of the second identity in Eq. 
(fTUl) ] upon exchanging z and z', 


J 2 =4 


dzdz' - 


rA[a2 In |Li(a2) + y - A + a^j - In \Li{p‘^) + A' - A + 


(15) 


and 
J 3 = 4 


f f dzdz'- 2 ^A[/3^ In 1,5^I - In |a^I]. (16) 

J — aJ ^ ^ 


Jq can be evaluated to be 

Jo = —8 In | 6 /a|[a 6 (l + 2 ln 2 ) In | 6 /a| + fc(3 + ln 2 )], (17) 


and accordingly 

I(k) = -— 81n|6/a|[a6(l + 21 n 2 ) In |5/a| 

16 k^ 

3 

+ fc(34-21n2)] + ^ Jj}. (18) 

i=l 


The Ji, J 2 , and J 3 are the same as the ones defined 
in Ref. [ 2 ^. The expression for Ji here appears to be 
different to that in Ref. [IJ but is actually equal. There 
is an extra term of the following form, 

/ dz [\z'^{X'-X)[L,{a^)-L,{f3% (19) 

J—a J ^ ^ 

in the expression for Ji in Eq. (10) of Ref. [ 2 ^, but 
this extra term can be shown to be equal to zero. (We 
note that in our foregoing derivation this term did not 
arise at all.) In fact, the factor of A — A' changes sign 
but that of Ti(a^) — Ti(/3^) [actually each of Li(a^) and 
Li{0^)\ remains unchanged upon exchanging z and z' as 
mentioned above. In consequence, 

/ dz / dz'^(A'-A)[Li(a2)-Li(/32)] = 0 . ( 20 ) 
J—a J ^ ^ 


Understandably, this term with null contribution might 
however cost a huge amount of labor in any further nu¬ 
merical or analytical calculations if not appreciated. 

The expression for I{k) has as a consequence been 
successfuly reduced in terms of two-dimensional integral. 
Similar type of reduction had been accomplished or ap¬ 
plied in Refs. [H, [2ll - l^ . [ 2 ^. (Brosens et al. [l^ in 
fact considered both of the static and dynamic cases.) A 
different approach was adopted by Holas et al. [l^ who 
essentially first calculated the imaginary part of Hi (k, w), 


and then the real part of it [with the static ni(k, 0) as 
its special case] via the dispersion relation, which virtu¬ 
ally also amounts to a two-dimensional integral over all. 
The remaining two-dimensional inte gral was then calcu¬ 
lated numerically in Refs. [^,[l^, ; analytically 
with the aid of some ingenious method and mathemat¬ 
ical assumption in Ref. l23|: and analytically but with 
computer programs in Ref. . We need point out that 
there exists discrepance between our Eq. m and the 
corresponding expression in Ref. [ 2 ^ [the one shown in 
the middle of the paragraph above Eq. (10) in Ref. [^]. 


III. EVALUATION FOR Ji 

It turns out that the integral in Eq. m for Ji can be 
fully carried through. We are about to achieve this step 
by step. Indeed, either the integral over z or z' can be 
carried out first. But before doing that it is worth point¬ 
ing out a somewhat hidden but useful truth. Ti(/3^) can 
actually be, by direct substitution in Eq. (nni), simplified 
to the following form: 

Li(/32) = 2|/3|. ( 21 ) 

Accordingly Ji, with some further algebra, can be writ¬ 
ten as 

Ji = 4 / dz^ f dz'^[a^R{z,z') + /3|/3|], (22) 

J-a ^ J—a ^ 


where 


i?(z, z') = Co(z)z'2 -f Bo{z)z' + Ao(z), (23) 


with 

Ao(z) = z^, Bo{z) = (2 -I- 2kz — k'^)z, C'o(z) = 1 - 1 - 2kz. 

(24) 

Eor brevity, Ao(z), Bo{z), and Co{z) will be denoted, 
respectively, as Ag, Bq, and Cq instead, i.e., with the 
argument z suppressed. We next define 

= [ dz—^a^jR{z,z'), (25) 

J — n. ^ 
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and 

Jf (z) = / dz'^m, (26) 

J —a ^ 

SO that we can write Ji as 



The integration over z' on the right hand side of Eq. (Ei) 
for Ji{z) presents no analytical difficulty. It is straight¬ 
forward and the resulting expression is, 

Ji{z) = 2kz — 1 — fc^/4 — z^(3 — 2 In \z\ + In |a 6 |). (28) 

On the other hand, the integration in Ji{z) is more sub¬ 
tle, but becomes also routine if we rewrite it as 

J(^(z) = z f dz'—^^jR{z, z') -|- f dz'\/R{z, z'). (29) 

J—a ^ J—a 

Each of the integrations in Eq. (Ei) can be carried 
through, and the result for J^{z) is 

J^iz) = 


+ 

where 


1 -I- -I- ^fcz -I- 2z^ -I- z^(2z -I- k) 

4 2 4Co 


zlzlIn 


aV'i(z) 


bip2{z) 


zBo 


A ^0 


•\/Co 


In 


i’siz) 


i’iiz) 


,(30) 


'ipi{z) — 2 Aq + Bob + 2\/A(fR[z, b), (31) 

tlj 2 iz) = 2Ao - Boa + 2y/A qR{z,- a), (32) 


^3{z) = 2^CoR{z, b) + 2Cob + Bo, (33) 

and 

Mz:)=2^/Q^R(z^-2Coa + Bo. (34) 

R{z, b) and i?(z, —a) can be written explicitly in the fol¬ 
lowing forms, 

i?(z, b) = [(1 -I- k)z + 6]^; i?(z, -a) = [{k - l)z -|- af . 

(35) 


Equation (l27l) , together with the explicit expressions (l28l) 
for Jj®(z) and ((5D1) for Ji{z), indicates that we have al¬ 
ready reduced Ji in terms of one quadrature. Further an¬ 
alytical refinement turns out to be practicable, but only 
after one finds a way to simplify the somewhat unwieldy 
quantities 'ipi{z)/'ip 2 iz) and ips^z)/'ip 4 (z). Indeed, we find 
the following simple expressions for them: 


tp2{z) 



2 

0 (z) + 



2 

0{-z), 


(36) 


where 9{z) = 1 for z > 0 and d{z) = 0 for z < 0; and 


i’sjz) 

tpiiz) 


V^-1 


(37) 


The simplifications play a critical role in our further re¬ 
ducing Ji to a fully explicit form [Eq. below]. The 
verification for Eq. Eni and Eq. (1571) is given in Ap¬ 
pendix A. We here substitute them into Eq. EQl) to ob¬ 
tain 


Jf^(z) = 1-1- -I- -kz + ( 2 — In 

4 2 ' 


1 - 


fc 2 


Bp 
4Co 
1 


(2z -|- k) 


+ -tttZ^P -k^ + 4fcz(6 - k^) + 12/c^z1r(z), 

(38) 


where we have adopted the following notation, 


Y (z) = In 


-I- 1 
y/C'p-l 


(39) 


We then substitute Eqs. (1751) and (1551) into Eq. (1771) and, 
with some further routine calculation, obtain 


Ji 


where 


20 + 161nfc -I- ^(fc^ — 1)^ In 
k 

{k^ - 2Ak + 8)ln\2b\ 

{k^ - 2Ak-8) ln|2a| 

i[37?i -H 2(3 - A:^)?7 _i - (fc^ 


fc + 1 
fc- 1 


-1)V3],(40) 


Vn = k f dzCl'^Y^z). (41) 

J —a 

The explicit expressions for 771 , 77 - 1 , and 77-3 are given 
in Eq. (|B2I) . Eq. (IB3|) . and Eq. (IB4I) . respectively. By 
substituting Eqs. (IB2I) . (IB3I) . (IB4I) into Eq. (1401) . we 
obtain our final result for J \: 

Ji = —16-I-4(5fc-I-4/A:) In | 6 /a|. (42) 

The result is surprisingly elegant. 


IV. EVALUATION FOR .h + J 3 

Among the three expressions as shown in (1751) . (I15L 
and (ITOl) . the expression in (1751) for Ji is likely the one 
the most amenable to analytical computations. This is 
in fact also the cause that it can be fully reduced to the 
explicit expression of (|471) . On the other hand, due to 
the logarithm forms in them, Eq. m and Eq. (HU) 
are much more difficult. Further reduction of J 2 and 
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Ja, at least that of J 2 , to one quadrature, had been 
believed beyond barehanded effort [l^. We find indeed 
that careless choices of procedures might easily make the 
task intractable. It is critical to search for optimal ones. 


Instead of evaluating J 2 and J 3 separately, we choose 
to combine them together and deal with them at one 
stroke. We write ^ 23 = J 2 +J 3 , and accordingly have 


•^23 = 4 / f dzdz' ^ In |[Li(q!^) + A' - A + Q!^]/a^| - /3^ In |[a - 6 + 2z + 2fi/\fi\]/(!i\) . (43) 

J-aJ ^ ^ 

We have also made the use of Eq. (EH) in obtaining the above equation. We next perform some manipulation on Eq. 
(l43l) to bring it into the following form: 

J 23 = 16 [ dz-X\n\4X\ [ dz'--m, (44) 

J-a Z z' 


where 


N = 


dzdz 'In |a^ + A' — A 


2./r(^)\ - ^2 In |/3(a -b + 2z- 2/3/|/3|)|]. (45) 

I- 


The first term on the right hand side of Eq. has vir¬ 
tually been reduced to a one-dimensional integral for the 
integration dz'\jz' is trivial. Hence we need concern 
ourselves only with the term of —4A^. For this purpose, 
we first rewrite N as the following: 


N = \ dz\\Ny{z)-N2{z% 


(46) 


where 
iVi(z) = 
and 


f<b 2 

dz'^ In la^ + A' - A - 2^/R{z,z')\, (47) 


N 2 {z) = j dz'^h\\2j5{z-b)\+[ dz'^\D.\2j5{z+a)\. 
J-a Z Z 

(48) 

In the preceding expression for N 2 , the integral over z' 
has been purposely separated into two regimes of —a < 
z' < z and z < z' < b. This is intended to take care 
of the subtle singular behavior of the logarithm factor 
In |/3(a — 6 -I- 2z — 2/3/|/3|)| in Eq. (|4^ . After performing 
partial integration on each of the integrals on the right 
hand sides of Eq. (H7)) and Eq. (H51) . one obtains 


and 


N 2 iz) = 2 I - — - zln 
ab 


In |2A| -I- 2zWi{z) In \z\ 


dz'I z' -^ — 2z In |z'| ] —. (50) 


Here we have introduced the following notation, 

z + a 


Wi{z) =ln 


z-b 


(51) 


Equation (l46)l . together with Eqs. (|4^ and ([50|) . is then 
substituted into Eq. (H41) to get 

J 23 = 8(2 In 2 — l)(a61n | 6 /a| -I- fc) In |5/a| 

-h4(-fc$i -h 24>2 - 2k<i>3 + 24>4), (52) 


where 


dzU'dz'^ ^ 


' —a ^ J —a 


z' yjR{z,z'y 


/ b A 7^ z' 

dz- / dz'—ln|z'|, 
.a z a/3 


(53) 


(54) 


iVi(z) = 2 1-- -I- zln 

' ab 


ln|2A| 


-b 


dz'l z' - - + 2 zln \z'\ 

z' 


1 

X — 


kz 


a\^R{z,z') 


- 1 


(49) 


and 


4)3 = / dzX [ dz' - . InU'l, (55) 

J-a J-a a^R{z,z') ' ^ ^ 


/ b \ 

dz-VEi(z)ln|z|. 

-a ^ 


(56) 
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$4 is clearly not our concern temporarily, since it is al¬ 
ready in terms of one-dimensional integral. Fortunately 
it turns out that the two-dimensional integral in can 
be fully carried through. We first attack it. To this end, 
we write 


$1 



dz-^i{z), 

z 


(57) 


where 


$i(z) = / dz' 


\/R{z,z'' 


■. — z 


dz' 


z'^/R{z,z') 


• (58) 


Each of the integrals over z' in Eq. (l58t is routine. Af¬ 
ter carrying out them, we obtain the following result for 


^i{z) = 




In 


aV’i(z) 


'<p3{z) 

bip2{z) 

•\/Co 

tp4{z) 


(59) 


or, with the aid of Eqs. dSi) and (EZl), 


l>i(z) 



y(z) -f In 


4a & 

”fc 2 " 


(60) 


Substitution of Eq. (| 6 Q[) into Eq. dig, followed by fur¬ 
ther algebra, yields 

< 1)1 = A:(ln |4a6| -2\nk) - + ^(1 + (61) 

The use of Eqs. (IB2I) and (IB3I) in the above equation 
then yields 

$1 = ^ [(8fc3 + 9 fc 2 + 4) In 6 -p (8fc3 - - 4) In |a| 

OK 

-ik + 16k^\n\2/k\]. (62) 


This is another elegant result indeed. 

^2 in Eq. (l54l) and < 1)3 in Eq. (1^ are our next con¬ 
cern. Clearly one can do nothing further directly with 
regard to integral over z' in those expressions, due to the 
simultaneous appearance of the factor In | z' | and nontriv¬ 
ial denominator in each of the integrands. Our trick is to 
reverse the order of integrations over z and over z'. This 
means to perform the integration over z first instead of 
that over z'. This consideration yields the following form 
for ^ 2 , 


< 1)2 = 2 / dzzln |z| 



z ap 


(63) 


Notice that we have exchanged the symbols of z and z' 
to obtain the preceding equation. The integration over 
z' can now be readily carried out. After that one obtains 


where we have introduced the notation 

z — a 


W 2 {z) = In 


z + b 


(65) 


This temporarily finishes our job for < 1 ) 2 , for it is already 
in terms of one quadrature. 

Next we consider < 1)3 of Eq. (1551) in the same fashion. 
First notice the property 


i?(z, z') = i?(z', z), 
which enables us to rewrite < 1)3 as 


$0 = 


dz<^z{z) ln|z|. 


where 


rb \/ 

$ 3 (^) = dz' 


ay^R{z,z') 


( 66 ) 


(67) 


( 68 ) 


Explicitly, 

$ 3 (z) = ( 6 -|-z)(a —z) / dz' 
+ {k + z) ( dz' 

J —a 


1 


-a ay/R{z,z') 

b ^ rb 


\/R{z,z') 


' dz' , 

\/R{z,z') 


(69) 


Each of the integrals in Eq. (1691) is routine. They can be 
carried through and the result is 


^3(^) = 


where 


A-[k + 2zf 
4fc|z| 


In 


<(' 1 ( 2 :) 


4>2{z) 


-W2{Z) 


+ \{k + 2 z )(2 -h 3kz)C~^^^ In 




Mz) 


2z “h k 
, (70) 


4>i{z) = 2k'^z‘^ — kz{2z+k){z—a)+2k\z\\/R{z, —a), (71) 
4>2iz) = 2k^z^ — kz{2z + k){hz) + 2fc|z| \/R{z, b), (72) 


(j) 3 {z) = 2\JCoR{z, b) + 2(6 -I- z)C'o - kz(2z -|- k), (73) 
and 

<() 4 (z) = 2y/CoR{z, —a) -|- 2(z — a)Co — kz{2z + k). (74) 

It is not difficult to verify that <(> 3 (z) = tjjaiz) and (j) 4 :{z) = 
V' 4 (z). Therefore, according to Eq. (1571) . 


4’3iz) f v^Cq -f 1 


(( 14(2 


VC~o-lJ ■ 


(75) 


r 1 

^ 2 = dz-ln|z| [2a61n | 6 /a|-I-AVFi(z) 

J — a ^ 

-{z^ + kz - ab)W 2 {z)l (64) 


On the other hand, it is shown in Appendix A that 


01 (^) 

02(2:) 



(76) 
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We next substitute Eqs. 123 and (l76)) into Eq. (l70)) , and 
then the resultant equation further into Eq. (|H7|. The 
result for $3 is in this way obtained as 


$3 



ln|z| 


2z + k 


{a 


z){b + z) 
kz 


W 2 (z) 


+(2z + k)(2 + 3kz)C~^''^Y(z) 


(77) 


By now we have virtually completed our task of reduc¬ 
ing J23 in terms of one quadrature. By substituting Eqs. 
(1551) . (1^ . (IMl) , and (1771) into Eq. (1521) . we reorganize 
our result for J 23 in the following form: 

J23=8[Po + (fc^-l)^"l + 2P2-fcP3], (78) 


where 


and 


P3 



dz(k + 2z){2 + 3kz)CQ ^^^Y(z) In \z\. 


(82) 


It turns out that the integral in the expression for P 3 can 
be further refined. Indeed, it is shown in Appendix B 
that 

P 3 = —(2fcln \k/2\ — 1/k) + 4(6In^ b — aln^ |a|) 

o 

-I- -(31n |2/A:| — 5 -|- 1/A:^)(61n6 — aln |o|) 

o 

“ 1)^1- (^3) 


Po = 


8fc2/31n|fc/2| -4/3 
[21n2 — I -|- In |a6|]a61n^ 1^/®! 

2[fcln2 — 6(4fc^ — 2fc -|- l)/3fc] In 6 
2[fcln2 — a(4fc^ -|- 2fc -|- l)/3fc] ln|a|, (79) 


Pi = 




(80) 


Remarkably, it can be seen from Eq. (1831) that the com¬ 
bination form of (fc^ — l)Pi — kP^, the form in which Pi 
and P 3 enter the right hand side of Eq. ([78]), can be fully 
reduced to an explicit form with no more integral in it. In 
this light we substitute Eqs. (IT^ and (1831) into Eq. (178)) 
to obtain the final expression for J 23 as the following: 


^23 = 8 b(fc) + 2P2], (84) 


r 1 

P 2 = / dz-[XWi(z) — (h +z)(z —a)W 2 {z)\\rL\z\, (81) where 

J — a ^ 

I 


g{k) = 2 -|- 4fc^ In \2/k\ + a61n |a 6 | In^ | 6 /o| -I- a6(21n2 — 1) In^ | 6 /o| -I- 4fc(aln^ |a| — 61n^ 6 ) 

-|-[2/c^ -|- 9fc/2 — 2/k — 2k(k -I- 1) In 2 -|- 4fc61nA:] In 6 

-|-[2fc^ — 9fc/2 -|- 2/fc -|- 2fc(l — fc) In 2 — 4fcaln fc] In |a|. (85) 


V. I(K) EXPRESSED IN TERMS OF ONE 
QUADRATURE; AND FURTHER IN TERMS OF 
SERIES 


Further substitution of Eq. (1421) and Eq. (|84|) into Eq. 
(fT 8 |) then leads to 

I{k) = -^[go{k) + 2 P 2 ], ( 86 ) 

where 

go(k) = 2k(a\T? |o| — 61n^ 6 ) - 1 - a61n |a 6 | In^ | 6 /a| 

— 2(61n6 — aln |a| — fc — A:ln |A:/ 2 |)^ 

-f 2k^(l + \n^\k/2\). (87) 


refined in the next section into the following final form: 




6 In" 


-a In 
3 


1 

X 

(1 ,, 

X 

,M 

+ / ax - T In 


1 — a61n 


— k 1 

Ja2 X-a^ 


V2 


J\ 


,(88) 


or, 


i(k) = 


2k^ 


a In^ 


-6 In 


-f a 


VP 


dx- 




\ 1 

1 —ab in 

V 2 

b^x 

— k jinx 


(89) 


As a consequence, we have accomplished reducing I(k) The resulting expression for I(k) is amazingly simple, in 
in terms of one quadrature. Equation (1861) will be further construct to its original six-fold integral form. Equation 





























Q in the Introduction is obtained from Eq. (1881) together 
with Eq. ([S]). 

The one quadrature in Eq. (18811 presents no numerical 
difhculty in applications. This is particularly so with the 

I 


aid of the expansion forms in the limiting cases obtained 
in Sec. VII below. But it turns out that we have even a 
better choice. In this section, a series representation for 
I{k) is developed. To this end, we rewrite Eq. (15^ as 


= 2 F 




h\n 

3 


abln 


In obtaining Eq. (19011 we have employed identities 

4. 

Jo 1 - a; 6 

and 

/•i 1 

/ dx- -ln^a; = 243). 

Jo 1 - a; 

oy series 
in Eq. (lOTl to obtain 


(91) 


(92) 


We now employ series representation 1/(1 — x) = 




fo(.k) + ^Uk) 


(93) 


where 


/o(fc) = aln^ 


a 

/i,, 

a 


( -6In 


b 

V3 

b 


+ a _-A: + a6C(3)-- 


6 


2b 


+ 


— ab + 2k In 

6 


2k 

1? 


In 



f a \ ^ 

,/l , 

b 

\ 


1+ 

k — ab\ —h In 



/ 

L \2b) J 

L V2 

a 

. 


(94) 


and 


fn{k) = a&ln 


+ k- 


ab 


1 


n + 2/ (n + 2)2 Vfe/ 

(95) 

We emphasize that the series on the right hand side of 
Eq. (IMl) converge for all fc > 0 and in effect the value 
of /(fc = 0) = — 1 can be equally (and rather easily in¬ 
deed) obtained from Eq. (1^ by employing the following 
identity, 


OO - 


(96) 


with C(2) = 7r^/6 explicitly. The series representation 
of Eq. (IM)l for /(fc) should be particularly useful for 
analytical purposes. 



0.5 


1 

k/2 


t5 


FIG. 1: The result for the quantity —7(fc), calculated with 
both of Eq. (15^ of this paper and Eq. (29) of Ref. [^ . 


VI. DERIVATION OF EQ. (USD AND PROOF OF 
THE EQUIVALENCE OF IT AND THE 
CORRESPONDING ONE OF EQ. (29) IN REF. 


It is by no means evident that /(fc) given as Eq. 
is equal to the corresponding one of Eq. (29) in Ref. 
( 2 ^. In fact the present author had certain reservations 
about Eq. (29) of Ref. [ 2 ^ for the mathematical as¬ 
sumption used in deriving it, and was convinced of its 
correctness only after a numerical check with Eq. (18611 . 
The numerical results from both of the expressions (be¬ 
fore the rigorous proof of their equivalence was found) 
are shown in Fig. 1. They attain the same illustration. 
Besides, it can be derived analytically from Eq. (I55D that 
/(fc = 2) = —7r2/24 [see also Eq. (11431) and Eq. (11441) 
below], which is the same as that derived from Eq. (29) 
of Ref. [ 2 ^. In this section we prove their equivalence 
by refining both of them into the form of Eq. 


We hrst refine the quantity given in Eq. (29) of Ref. 
[ 2 ^ into the form of Eq. (I55D . For the purpose of conve¬ 
nience, we call it I^^{Q) (with Q = kj2). To this end. 
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we make the use of the following identity, 


where 


^ 1-x^ 

dx -— 


In^ 

1 + X 

1 — X 


rQ 

+6 

/ dx 


Jo 

of Ref. 1 


= -(Q+Q)ln 


1 + Q 




l-Q 

1 + a; 


1 — a; 


, (97) 


fQ 

I^^{Q)=UJl +L 02 dx 

Jq 


1 + a;^ 
rQ 1 _ 

+ UJ 3 dx -^—In^ 


In" 


1 + a; 


1 — a: 


1 + a: 


1 — a: 


, (98) 


where 


1 — , 3 1 + Q 


l-Q 


, 0J2 = 


1-Q^ 

4Q2 


(99) 


sVg^ 2Q2 


1 + g 


i-g 


( 100 ) 


We next make the variable transform x = {y — a)/{y +a) 
in Eq. (|M|l to obtain 


I^'^{Q)—wi+ f dy [ uj 2 {y‘^ — a'^) + Sojsa^y^ ] 
J a 


1 


y{y^ - 


In" 


( 101 ) 


Further algebra then leads to 


I^^iQ) = ui- 
rb^ 


UJ2 , 

in 

b 

2 a^L 03 

In^ 

b 

[ 3 

a 

k 


a 


dy 


In 

y 

2uj3a^ 

^ In 

y 

4 “ 


1 

y \ 

y-a^ 



( 102 ) 


The remaining derivation to obtain the form of Eq. (IMll 
is routine. 

We next refine I{k) of Eq. (|Mll into Eq. (IMll . To this 
end, we carry out partial integration on the right hand 
side of Eq. m to rewrite it as 

P 2 = A:[61n^ 6 In \2b/k\ — aln^ |a| In |2a/fc|] + P 2 , (103) 


with 

P 2 =-z f dz[(2z — A:)W'i(z) + (2z + A:)W2(2;)] In^ |z|. 

4 J-a 

(104) 

With some straightforward algebra, P 2 can be put into 
the form 


P 2 = (g - m 2 - m 3 - t 4 )/ 8 , (105) 


j-b^ 

ti= dy\n^y\n\y - a^l, (106) 

Ja^ 


rb ^ 

1 dv In" V In 

y + a 

J\a\ 

y-a 


(107) 





dy In^ y In 


y + b 
y-b 


and 


U= dyln y\a.\y - b^\. 
Ja‘^ 


(108) 


(109) 


The next manipulation is the key to the procedure. For 
ti we perform further partial integration to obtain 


ti = [ 6 ^((i 4 + 1 ) - a^(di + 1 )] In | 2 /c| 

. 6 " 1 

- / dy -^[y(ln^y- 21 nj/ + 2 )-a 2 (d^ + l)], 

Ja^ y- 

( 110 ) 

where 


di = 21 n|o| — 1 , (i 2 = 2 (ln |a| — 1 ), 

da = 2(1 —ln 6 ), d 4 = — 21 n 6 +l, (HI) 


with d 2 , ds for later quotations. Further algebra then 
yields 

ti = 6 ^(d 4 +l) In | 2 A:| —(d 4 + l)( 2 A:+a^ In | 2 A:|)—Si, ( 112 ) 


where 


/ dy ( 2 d„ + In 

JL 

)ln 

]L 

/a2 y - a" V 


J 



(113) 


with S 2 , ^ 3 , and S 4 for later quotations. Similarly we 
can obtain 

g = ^[b{dl + 4) In \2/k\ — a{d\ + 4) In |2a^/fc| + aS 2 ]. 

(114) 

For ti we first perform partial integration to get 

ti = [b^{dl + l)-a^{dl + l)]\n\ 2 k\ 

1 

— / dx --r[a:(ln^ X — 21 na: + 2 ) — ^^(d? + 1 ) 1 . 

Ja 2 x-¥ 

(115) 

We then make further variable tranform x = a^b^/y to 
get 

ti = (d 4 + 1 )( 6 ^ In \ 2 k\ — 2 k) — a^{df + 1 ) In \ 2 k\ + a'^b'^Ei. 

(116) 
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Similarly we can obtain 


where 


ts = ^a{d\ + 4) In |fc/2| + 6(^3 + 4) In \2b'^/k\ — 0 ^ 623 ], 

(117) 

Substitution of Eqs. (11121) . (11141) . (I116L (I117p into Eq. 
(I103|) yields 


r = [2i + ka'E.2 - kba^Es + a'^b'^Eil/A, 


P 2 = fc(ln 6 -l)(ln5 + 261n|25/fc|) 

— fc(ln |a| — l)(ln |a| + 2 a In | 2 a/A:|) — T/ 2 , 

(118) which has the following explicit expression: 


(119) 


J 


1 

T = [2A:ln \ab\ + b{b — 2a + 2a In 6 ) In |t>/a| — 2k — 2a6/31n^ l^/o|] In \b/a\ + / dx -^ In 

J ^ ^ 

( 120 ) 

One then substitutes Eq. (I118|) together with Eq. (11201) into Eq. (IM)) for I{k). The remaining derivation to get the 
form of Eq. (1551) becomes routine. 



VII. EXPANSIONS IN THE LIMITING CASES 


The limiting structures of ni(fc, 0) at fc = 0, fc = 2, and at large k have various physical implications and deserve 
particular attention. In virtue of our proof of the correctness of the main conclusion of Engel and Vosko’s theory in 
Ref. 1^. th ere is no doubt that the expansions obtained by them in the limiting cases [i.e., Eqs. (30), (31), and (32) 
in Ref. [^] must also be right. Indeed they are all confirmed with calculations based on Eq. (1551) . We sketch the 
calculations below. 

It turns out that the expression (|89D is rather suitable for deriving the expansions. We first consider the cases of 
fc —>■ 0 and fc —>■ 00 , To this end, we make the variable transform x = [(1 — y)/(I + y)]^ in Eq. (1551) and bring it into 
another form. 


I{k) = 


2 /c 2 


6 In 


: In 


2 Pdyf- 

Jo \y 


1-2/ 


In 


1 + y 


1 - 2 / 


ab In 


a(I + y) 


K^-y) 


- k 


We act the operator of ^(1 — Q^)-^q + 1 on the quantity I{k)Q‘^. With some algebra, one can obtain 


( 121 ) 




[mQ^] = 


4Q I 



b 

, fi 2 \ 

1 + 2 / 

blxF 


- 2 / dy[- + - -) In 


. 

a 

io V 2 / 1 - yJ 

1 - 2 / J 


( 122 ) 


One then acts further the operator of {1 /Q)d/dQ (which 
is of course just 2d/dQ'^) on the both sides of the preced¬ 
ing expression. This action yields the following result. 


d 

QdQ 




d 

QdQ 


-bl 


[I{k)Q^] 


(123) 


where 


F{Q) 



+ 





In 


1 + y 
1 - 2 / 


• (124) 


It can be readily shown that 


F{Q^0)=4Y. 


1 


n—0 


and 

00 

F{Q 00 ) = 4^ 
where 


n—0 


2n -|- 3 


2n + 3 


24'(n) 


1 


2n + 3 


- 1 


2^'(n) -I- 


1 


2ti F 3 


- 1 


4'(n) = 


1 


i=0 


21 + 1 


Q2n+3^ 

(125) 
1 

Q27T3> 

(126) 

(127) 





































































In passing we note the remarkable symmetry in the above 
two forms. For fc —>■ 0 we have, from Eq. (I123p and Eq. 

(insi), 


QdQ 

oo 

= y— 

0.71 -\ 






n =0 


2n + 3 


2 «'(n) 


1 


2n + 3 


- 1 


Q^’". (128) 


With Eq. (11281) . it is straightforward to determine the 
expansion for /(/c —^ 0). To this end, we write 


/(A:^0) = -l + ^c„Q"", (129) 

n— 1 

where c„ are the coefficients we have to determine. To 
write the preceding form, we have made the use of the 
truth of I{k = 0 ) = —1 which can be readily obtained 
from Eq. (11211) and is in fact well known [ll|, [2]] - l^ . The 
fact that the expansion retains only even powers of Q is 
evident from Eq. (11281) . [In fact I{k), which is defined 
physically only for A: > 0, is an even function if extended 
to the range of A: < 0. This fact can be best appreciated 
from Eq. (1551) . but also in another interesting way. Upon 
making the variable tranform x = 1/y in Eq. (j89ll . one 
can bring it into yet another form, 


1 


h 

fi 

b 

\ 


b\n^ 

a 

(r'” 

a 




a^x 

A , 1 



— k \ in X 


(130) 


But, on the other hand, the right hand side of the above 
equation can be exactly obtained from that of Eq. (|M1) 
by merely changing k to —k. Thus I{k) = I{—k).] We 
then substitute Eq. (11291) into the left hand side of Eq. 
(I128|) . By equating the coefficients of the same powers 
of on both sides of the resulting equation, we obtain, 
for n > 1 , 


(n + l)(n + 2 )c„+i — n{n + l)c„ 


2(2n + 3) 


2 «'(n) 


1 


2n + 3 


- 1 


(131) 


or further 


n(n+l)c„ = 




1 


2^21 + 1 
1=1 


24'(Z-1)- 


1 


2 / + 1 


■-1 


■ (132) 


By the use of Eq. (|127p . we have 


1 

n(n + l)c„ = - - l)]['k (0 + d/(/ - 1 ) - 1 ]. 

( 133 ) 
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After carrying out the above summation we have, with 
the aid of the fact 'I'(O) = 1, 


'I'(n)['k(n) — 1 ] 
2 n(n + 1 ) 


(134) 


and accordingly 


I{k 0 ) 


^ , l^ vl/(n)[4^(n)-l] ^,„ 
2 ^ n{n + l) ^ 

n—1 ^ ' 


(135) 


Equation (jl35l) is the same as Eq. (30) of Ref. 

In particular, the leading two terms in Eq. (jl35l) which 
are 


I{k —>■ 0) — —1 + gQ^ T (136) 

had been obtained numerically by Kleinman and Lee [2l| 
and later by Chevary and Vosko (who also obtained 
the next order 46(3^/675 correctly). The form in Eg . 
(11361) y ields the following leading gradient correction [la . 
I28l432l | to the local densit y appro ximation of Kohn-Sham 
exchange energy density [ 2 ll - [^ 


A 41 ) = 


5e^ 

2167r(37r^)iA 


[V»(r)] 

[n(r)]‘*A 


2 


(137) 


The leading gradient correction to the Kohn-Sham ex¬ 
change energy was investigated originally by Sham [s^, 
but he got the coefficient to be —7e^/[4327r(37r^)^/^] in¬ 
stead by employing screened Coulomb potential with 
screening going to zero (see also Ref. [s^). Later Gross 
and Dreizler [33 (see also Ref. (s^) got the same result 
as Sham’s, essentially employing also screened Coulomb 
potential (with screening going to zero). This well-known 
controversy has been one of the major causes inspiring 
the investigations of the exact structures of Hi ( k, 0), and 
has been effectively elucidated in Refs. [2ll - l^ . Our re¬ 
sult based on the rigorous derivation helps to give a final 
settlement of the controversy itself. Questions such as 
whether in Eq. (11371) or the one by Sham should 

be added to the type of Ma-Brueckner’s gradient cor¬ 
rection to the correlation energy (to get the Kohn-Sham 
exchange-correlation energy be yon d the local density ap¬ 
proximation) remain open [2ll l31j . 

For A: —>■ 00 , we have, from Eq. (11231) and Eq. (I126L 


d 


QdQ [2 


= E 




1 


n=0 


2n + 3 


QdQ 
24'(n) 4 


+ 1 


[Iik)Q^ 


1 


2n + 3 


- 1 


Q2n+e 


.(138) 


It is not difficult to see from Eq. (|138l) that the expan¬ 
sion of I{k) for large k commences with 0 ( 1 /Q®), and 
accordingly it assumes the following general form. 


OO 

J(fc -» 00 ) = c„ 

n—3 


1 


( 139 ) 
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Substituting this form into Eq. (jl38L we obtain 


(n + 2)(n + 3 )c „+3 - (n + l)(n + 2 )c „+2 


2(2n + 3) 


1 - 


1 


27z -f 3 


- 2 «'(n) 


(140) 


for n > 1. Similar procedure to that used above for the 
case of A: —>■ 0 can be employed to obtain 


4'(n-2)[^'(n-2)- 1 ] 
2(n — l)n 


(141) 


Thus 


I(k ^ oo) 


2 ^ (n + l)(n + 2 ) Q 2 ("+ 2 )’ ^ ’ 

n—1 ' ^ / -v 


confirming Eq. (31) in Ref. [^ . (there is a sign error 
in the first identity of that equation.) We note that the 


leading two terms of Eq. (11421) had been reported by 
Geldart and Taylor dl (see also Ref. 0 ) ; the leading 
term had also been obtained by Kleinman [2^, . 

We next give the expansion near k = 2. With the aid 
of Identity (|TO . the characteristic property of /(fc = 2) = 
—7r^/24 and the leading singular term (meaning the most 
divergent term upon derivative with respect to k) can be 
immediately seen from Eq. (I89p to be: 


I{k 2 ) = 




a 

2 


(143) 


The conclusion drawn in Ref. that the singularity of 
ni(k, 0) dominates that of no(k, 0) is corroborated. In 
general, we obtain 


J 


I{k ^ 2 ) = 


If 7r2 


+ ( 2C(3) —:r "*■ ( 3C(3) —+ o 1®^ + ( 4C(3) “ 777 + ( 5C(3) + 


6 2 


24 


Il7r2 37 


144 96 


+ 6C(3) + 


o'" 


377r2 125 Y 

192 ~ T^J 

- 00 


n—1 


00 

^ e„a" In^ 

n—0 


- 1 )a + ^ Tn+lO^ 

n=2 


In 


, (144) 


where 


= 2 (n+l)+^ 


1=0 


I — n 1 
{l + l){l + 2) 2 ^+ 1 ’ 


(145) 


and 


Tn = 


-(n+l)+7(n-l)-80(,.-2)+^ (i"-^ 


■)3—m 


m—3 




3(m-1)2-1 ™ 1 


li_V-2' 

i{m — l)(m — 2 ) ^ 


, (146) 


with (j){m) = X)™ol/(^ + 1 ) ■ particular, eo = 2 , 
Cl = 15/4, 62 = 131/24, and 63 = 229/32; r 3 = 277^/3, 
r 4 =57r2/6+13/8, and =7r2+173/48. Equation (11441) 
confirms Eq. (32) in Ref. [ 2 ^ in virtue of the following 
identities. 



= 3C(3), (147) 

where C(3) = 1.202056903.... We would like to point out 
in this connection that Eg . (21) in Ref. is at variance 
with Eq. (32) in Ref. and accordingly also with Eq. 


(11441) above. 

Appendix A: Verification for Eqs. (1361) . (1371) . and 

dZU 

It is worth pointing ont first the following fact, 

(1 + fc)z + 6 > 0 , and (fc—l)z + a> 0 , (Al) 

which can be readily verified for —a<z<b. Accord¬ 
ingly, 

\/R{z,b) = {l+k)z+b; ^/R{z,—a) = {k—l)z+a. (A2) 
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A direct substitution of Eq. (|24)) and Eq. (IA2I) into Eq. 
m and Eq. (15^ . respectively, leads to 

'ipi{z) = Ab'^z{z + a), and = k'^z{z + a), (A3) 

for z > 0 ; and 

V'i(z) = — 6), and ip 2 iz) = Aa‘^z{z — b), (A4) 

for z < 0. Equation (IMll is hence verified. Equation (|76ll 
can be verified in the same manner. 

We next verify Eq. dSZl). Once again we make a direct 
substitution of Eq. (l24t and Eq. (IA2I) into Eq. (l33)) and 
Eq. dMl), respectively. By carrying out some algebra, we 
can obtain 

tpaiz) ~ [ C'q + 2 {k + 1)00^ + k{k + 4)Co 


and 


+ 2{k^+ k-l)y^ + k^- 1 ], (A5) 


and 


1 


= [ C^ + 2{k-l)C^o^^ + k{k-4)Co 

+ 2{-k^ + k + l)^/^ + k^-^]. (A6) 

Both of these forms can be factorized, and the result is 

Mz) = ^(v^+l)'(v^+fc+l)(v^+fc-l), (A7) 


Mz) = ^(v^-l)"(v^+fc + l)(v^+fc-l). (AS) 


Equation dSzl) is as a consequence established. 


Appendix B: Evaluation for rjn in Eq. HUD and 
derivation for Eq. 


An efficient way to evaluate r]n and P 3 is to introduce 
the variable transform, z = ^(y^ ~ l)j which casts Eq. 
(I 4 T]) in the form 

l-k 

Vn= / dy(y + l)"+iln 


y + 2 


(Bl) 


The remaining calculation for rjn is then routine. The 
explicit expressions for 771 , rj-i, and 77-3 are given as 
follows. 


rji = —[2k + {k^ + fc + l)61n \2b\ — (k^ — k -\- l)aln |2a| — k{k^ + 3) InA:], 
3 

77_i = 2(6In |26| — a In |2a| — /c In k), 


and 


??-3 = 2 


— 1 


In fc + In 


k + \ 


k-l 


In |26| + --- In |2a| 


k + 1 1 — k 

1 /'„,2 


We next evaluate P 3 to verify Eq. (IM|l . With the transform z = ^(y^ — 1), Eq. (1^ can be rewritten as 


2fc2 Z! 


nl+fe 


In 


dyy'^^ In 


y2 - 1 


2 k 


In 


y + 1 


y- 1 


n=-l •7|l-fc| 

where 

7i =3, 7o = - 2, 7 _i = - 1. 

One then carries out partial integration to bring Eq. (IB5|) into the following form, 

1 


P 3 = 


2fc2 ^ 2n + 1 

n— — l 

with k = |1 — A:| and defined as 


(l + /c)^"+On61n 


26 

-P"+Mn|a|ln 

fc + 1 


k 


k-l 



p\-\-k 

^n= dyy^”’^^ 

Jk 


.y +1 


■In 


2 k 


(y +1)^ 


y- 1 


■In 


2 k 


{y-if 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 
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for n = 1, 0,-1. The ll-i can be evaluated as 


Pb ^ 

f2_i = ho{k) — 2k dz— In \z\ 

J-a ^0 


where 


For n = 0,1, one has 


or 


ho{k) = In \2k\ In 


+ In^ \2b\ - In^ + 1| + In^ k - In^ \k - 1|. 


^r?. — 


p2b 

/ dx{x — 1) 

Jk+l 


2n+l 


■ In 


2k 

™2 


[ dx{x + l)^'^+^-\n 

Jk-i X 


2k 


2n+l r p 

^n=Yl C2n+1 i-r I 

_1 . d k 


2b 

fe+1 


dxx'^ ^ In 


2k 

t ,2 


dxx^-^ In 


/fc-i 


2k 

T2 


By carrying out the integrals on the right hand side of Eq. (IB12I) . we obtain, for n = 0,1, 


2n+l r / ^,2 \ / h \ 

fl„ = /io(fc)+ ^ C 2 ™+i^(-r (25rfmln— - 2 j + (-fc)™ f m in - - 2 j 

m = 1 L\ / \ / 


+(k + 1)™ ( m in 


(25)™ ( min 

2k 


{k + iy 


+ 2 j + (1 — fc) 


-2) +(-fc)™( min 

2k 


Tn 


{k-iy 


+ 2 


(B9) 


(BIO) 


(Bll) 


(B12) 


(B13) 


One then substitutes Eqs. (IB9I) and (IB13I) into Eq. (|B7|) to obtain, with some further algebra, Eq. (1551) . Or one can 
obtain first explicit expressions for Oq and Oi. In that case, one has 


^0 = ho{k) + 2[A:lnA: — 25In |25| + (fc + 1) In |fc + 1| — (fc — 1) In |(5 — 1|], 


(B14) 


and 


Oi = /lo(fc) + - [ —20fc — 12fcln2 + k{6 + 9fc + 2k^) Infc — 25(8 — k + 2k^) In |25| + (fc + 1)(11 — 5fc + 2fc^) In \k + 1| 


— (fc — 1)(11 + 5fc + 2fc^) In — 1| 


(B15) 


In any case, the remaining derivation to obtain Eq. (|83ll 
is straightforward. 
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